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Abstract 

Event rates of various lepton flavor violating processes in the minimal super- 
symmetric SU(5) model are calculated, using exact formulas which include Yukawa 
vertices of lepton-slepton-Higgsino. We find subtlety in evaluating event rates due to 
partial cancellation between diagrams. This cancellation typically reduces the event 
rates significantly, and the size of the reduction strongly depends on superparticle 
mass spectrum. 
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Lepton Flavor Violation (LFV) is an indication of physics beyond the standard model. 
Supersymmetric (SUSY) extension of the standard model, which has been strongly moti- 
vated as a solution of the gauge hierarchy problem, generally possesses LFV as generation 
mixing of slepton masses [TJ. In particular, interaction of superheavy particles may induce 
LFV in the slepton masses through renormalization effects, giving non-negligible event 
rates for LFV processes [[|. Concrete realization of the mechanism has been considered 
in supersymmetric grand unified theories (GUTs) such as 577(5) @, [3], f|, 5*0(10) (|, ^ || 
and flipped SU(5) 0, as well as in models with right-handed neutrinos || |9], liPf . 

A remarkable observation has been made in Refs. || [|, which has pointed out that 
large top-quark Yukawa coupling yields large rates of LFV processes in SUSY GUTs 
only one or two orders of magnitudes below present experimental bounds. The large top 
Yukawa coupling gives a sizable flavor violation in slepton masses squared through the 
renormalization effects above the GUT scale. Importantly || |J, the amount of the off- 
diagonal elements in the slepton mass matrix has mild dependence on the value of the 
top Yukawa coupling even when it closes to its quasi-infrared fixed point. 

However, this stable prediction of the LFV source in the slepton masses will not 
necessarily lead to a stable prediction of event rates. The purpose of this paper is to give 
exact event rates of fi — > ej, \i — > eee, fi-e conversion in nuclei, and r — > /i7 in the case of 
the minimal SUSY SU (5) model, and show that this is the case. When calculating these 



event rates, we include Yukawa type vertices of lepton-slepton-Higgsino [10|, which was 



not taken into account in Ref. f4|. As was already suggested in Ref. p0| , new diagrams 
coming from the Yukawa vertices interfere destructively to the other dominant diagrams, 
leading to reduction of the event rates. We find that this partial cancellation is subtle 
and generally has complicated dependence on superparticle mass spectrum. It takes place 
significantly under the universal scalar mass hypothesis because contribution from each 
diagram is comparable, and event rates we obtain are typically much smaller than the 
ones previously analyzed || |[]. Moreover, we will show that, given a LFV scalar mass 
matrix, variation of other parameters within a reasonable range will change the event 
rates by one order of magnitude or more. Without detailed knowledge on superparticle 
mass spectrum, it would be difficult to predict these LFV event rates. 

First of all, let us introduce the model we consider. In the minimal SUSY SU(5) model 
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both quarks and leptons are embedded in the 5* and 10 representations. The left-handed 
(SU(2) L doublet) leptons are in 5*'s, while the right-handed (SU(2) L singlet) ones are in 
10's. The Yukawa interaction in this model is described by a superpotential 

W = IffViVjH+VlffViQjH (1) 

with and fl being Yukawa couplings. Here supermultiplets ^ and correspond 
to the 10 and 5* representations, respectively, and indices % and j = 1,2,3) stand 
for a generation. H (H) is a Higgs supermultiplet in 5 (5*) representation. In the 
decomposition of the standard model, /„ becomes Yukawa coupling matrix of the up-type 
quark sector, whereas /a involves those of the down-type quark sector and of the lepton 
sector. The Yukawa couplings in Eq. (|J) break lepton-flavor number conservation as well 
as its quark counter part: one can diagonalize one of the Yukawa couplings matrices, fd for 
example, with respect to generation indices. When it is done, the other Yukawa couplings 
/ u among the 10's cannot be diagonalized in general and thus yield flavor violation. In 
the model we are considering, we assume that this is the only source of LFV. 

Above the GUT scale Mq k2x 10 16 GeV, soft SUSY breaking terms for scalar fields 
are given by 

- £ SU SYbreaking = ( m lo)^ t Vj + ( m l) i 4> jl 4> j + mffih + TH^h h 

+ (^A^ j h + V2A i ^ i cj )j h + h.c.), (2) 

where we have used lowercases ipi, (pi, h, and h to denote scalar components of the 
supermultiplets $j, H, and H, respectively. Soft masses for right-handed sleptons are 
identified with m\ and those for left-handed sleptons are m\. At the (reduced) Planck 
scale, M ~ 2 x 10 18 GeV, it is assumed that the matrices m\ and m\ are diagonal and 
A u and A d are proportional to the Yukawa couplings / u and /d, respectively, in order to 
avoid too large flavor changing neutral currents. When the energy scale goes down below 
the Planck scale, they will acquire flavor-violating components through renormalization 
effects which involve the flavor-violating Yukawa couplings. The magnitude of the LFV 
in the soft masses can be evaluated by using renormalization group equations above the 
GUT scale. At one-loop level, some of the equations are written as 
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+(10/ d /j + 3/ u /t)^ + (8A d f\ + 6A u fl)\ff. 
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(3) 



(4) 



(5) 



(6) 



In these equations, g 5 is the gauge coupling constant of the unified gauge group, M 5 is the 
corresponding gaugino mass and Q stands for the renormalization point (energy scale). 

Solutions of the one-loop renormalization group equations above the GUT scale were 
discussed already J4j. Following Ref. [f|], let us take the basis where the Yukawa matrix / u 
is diagonalized. Dropping the Yukawa couplings and A-parameters except for f^ 3 (= ft) 
and A^ (i.e. those for top quark) in the renormalization group equations, one can solve 
them analytically.^ It turns out that the universality for the flavor indices no longer holds 
in m 2 , A u , and A d , while it holds in m 2 . In particular, the matrix m 2 at the GUT scale 
is generally written as 



/ m 2 


V o 
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(7) 
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where m 2 is a generation blind contribution and / stands for a contribution from the top 

Yukawa coupling. Below the GUT scale, it is more convenient to use the basis in which 

the Yukawa matrix for leptons are diagonalized. Such a basis can be obtained by using 

the Kobayashi-Maskawa (KM) matrix at the GUT scale. In this basis, soft masses for 

*The renormalization effect by the bottom Yukawa coupling decreases flavor-diagonal, third generation, 
masses of the left- and right-handed sleptons, but does not significantly induce LFV masses. This could 
be easily understood in the basis that / d is diagonal. Though we will include the non-zero bottom Yukawa 
coupling in our numerical analysis, for the moment we ignore it for simplicity. 
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the right-handed sleptons exhibit LFV and their amounts are proportional to a bilinear 
of elements of the KM matrix (V): 



[m 



2 \j 

Ir> * 



m 2 S{ - VfVjl (Q ~ M G ) 



If universality of soft masses is imposed at the Planck scale, the quantity I is compactly 
expressed in the minimal SUSY SU (5) model as ||] 

/ = p G [m 2 + -{l-p G )A 2 + 0.198(1 - PG )A M 5G 

+ (0.224 - 0.029p G )M 5 2 G ]. (9) 

M$ G is a gaugino mass at the GUT scale. Precise meaning of an universal scalar mass ttiq 
and an universal A parameter Aq at the Planck scale will be given below. In the above 
equation, p G is defined as 

where f™ x (M G ), which equals 1.56 in the minimal model, is a maximal value of the top 
Yukawa coupling at the GUT scale under the assumption that the theory is perturbative 
below the Planck scale. Conversely if one starts with a large value of the top Yukawa 
coupling at the Planck scale, it converges to f™*{M G ) at the GUT scale and hence it 
exhibits a property as a (quasi) infrared fixed point. Thus p G < 1 and equality holds 
if the Yukawa coupling is on the fixed point. The structure of Eq. (|9]) is simple. It is 
proportional to p G as a whole, and the quantity in the parenthesis is of the order of 
the SUSY breaking scale. Furthermore as the top Yukawa coupling gets close to the 
maximal value, I converges almost to a value m 2 ,. This simplicity is an important finding 
ofRefs. [|,|. 

For the reader's convenience, we present in Fig. 1 dependence of the LFV source on 
the top quark mass. The horizontal line is the top pole mass (m]? ole ) divided by sin/3, 
which is defined as tan/3 = (h-i) / (h\) with h\ (h 2 ) being the Higgs boson that give 
masses to down-type quarks and charged leptons (up-type quarks). The vertical line 
represents (ml )e/mf where (m~ is a (p, e) component of the right-handed slepton 
masses squared, and m^ R is a physical right-handed selectron mass. These are evaluated 
at an electroweak scale, i.e., at the Z-boson mass. Here we choose M§ G and mo such that 
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a Bino mass M\ is 50GeV and mg fl = 300 GeV, and Aq = 0. The real line, dashed line, 
and dash-dotted line correspond to the case for tan/3 = 30, 10, 3, respectively. 



Complete formulas of the event rates for the LFV processes were given in Ref. flC 
where we took a mass-eigenstate basis for the sleptons. Instead of repeating them, we 
would like to discuss the qualitative behavior of the LFV rates using an approximate 
mass-insertion formula. In the approximation, the mass-eigenstate basis for the leptons is 
taken, and LFV masses squared as well as left-right mixing slepton masses are treated as 
perturbation. (This approach gives us a good approximation when the generation mixing 
in the chirality conserving scalar masses dominates over that in the A parameters. This is 
checked numerically in the model we are considering.) In this discussion, we concentrate 
on \i — > ej. The process is described by effective electromagnetic-dipole type matrix 
element: 

T = ee a *u e m^a al3 q^A^P L + A*P R )u^ (11) 

where Pr/l = (1 ±75)/2, is a muon mass, e the electric charge, q a photon momentum, 
and e a is the photon polarization vector. From this equation, the event rate of /x — ► ej is 
given as 

T(^e 1 ) = ^ml(\A^ + \A^). (12) 

l07T P 

In the model we are discussing, A^ is almost zero since left-handed sleptons have only 
tiny LFV soft-breaking masses with the non-vanishing bottom Yukawa coupling. The 
diagrams contributing to this matrix element have to have chirality flip of lepton and, to 
be proportional to SU{2) L x U{1) Y breaking vacuum expectation values. Such diagrams 
in this model are shown in Figs. 2. Fig. 2(a), in which the chirality flips on an external 
line, gives the following contribution to A\ ,[] 

lay 1 (m~ )% 
6 47r mf mf 

Ir Lr 

where, m,~ is an averaged right-handed charged slepton mass, (mf a (/i, e) component 

R iff 

of the right-handed charged slepton mass squared matrix, and ay = g Y /^ with gy 

^In our case, diagrams with two mass insertions like (m~ )^(w? ) T e could give a non-negligible con- 
tribution. However, we do not incorporate such higher terms since they do not change our argument 
drastically. Especially they do not change the sign of each contribution discussed here and subsequently. 
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being the U(l)y gauge coupling constant. To make our points clearer, here we adopt an 
expansion with respect to M\ /m 2 - , though we will show the exact results in the following 

''R 

figures. This expansion is motivated by the fact that the right-handed slepton masses 
become larger than the Bino mass through the renormalization effects, unless sfermion 
masses are negative at the Planck scale. (For example, m^. R ^ 2M\ in the case of the 
minimal SUSY SU(5) model.) 

In Figs. 2(b), (c) and (d), the chirality flips on an internal line. Fig. 2(b) gives the 
following contribution, 



i 2 

^2 lb = ^T-=2 =2 =2 • l i4 J 



la Y 1 (mjJZM^AJfp + ntmP) 
2 47r to- to- m~ 

Ir Ir '£ 

In this equation TOY is an averaged left-handed charged slepton mass, a SUSY-breaking 
trilinear coupling proportional to muon Yukawa coupling constant f^. It is notable that 
A 2 |b grows as tan/?, and can dominate over A 2 | a if 3Mi/itan/?/m~ ^ 1. Remember that 
matrix element Eq. ( p]) is proportional to the Yukawa coupling of the lepton, and to one of 
SU (2)l x [/ (l)y breaking vacuum expectation values, (hi) or (h 2 ). Some diagrams hitting 
the vacuum expectation value of h 2 result in the enhancement if tan j3 is large [|| [H] . 



Once the Yukawa couplings of lepton-slepton-Higgsino are correctly included, Figs. 2 
(c) and (d) also contribute to /i — > as 

a l I ay 1 K fl )g Afr/x tan/3 // 2 

47T mf mf mf mf 

Ir l R. IR l R 

A L \ aY 1 im V* f( » 2 ) (m) 

A 2\d ~ 73=2 =2 J2l=2-J, (J-Dj 



< 2 2 J « V 2 

47T mi 

(J? <ii Ir 



with 



o — 11a; + 4x 2 — x 3 + 2(2 + x) \ogx 

A(*) = ^rr^i > ( 17 ) 

1+ 4a; - 5x 2 + 2x(2 + x) log x 

T2\x) — ; h • (18) 

2(1 -xY 



Functions fi(x) and f'2{x) are monotonously decreasing functions of x and positive- 
definite. As x goes infinity, both functions vanish. Thus A 2 \ c and A 2 \d are not neg- 
ligible unless | /i | is much larger than to z ~ . As a result, A 2 |b and ^4^ | c are important when 
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Mi/^tan/3 is large, otherwise A\ | a and A\ |d will dominate. An important point is that the 
relative signs between A\ | b and | c , and that between A\ | a and Id are both opposite. 
This fact implies possible cancellation between the diagrams. We will soon show that this 
indeed occurs. 

We are now at a position to present our result of numerical computation using the 
formulas in Ref. |l(J. First, we consider the case where all scalars have a common soft mass 
at the reduced Planck scale (i.e., a universal scalar mass). The model is characterized by 
the following parameters at the scale M: a universal scalar mass tuq, a SU(5) gaugino 
mass M5, a universal trilinear scalar coupling Aq, a mixing mass parameter of the two 
Higgs scalars B, and a supersymmetric Higgsino mass parameter fi. Then, the boundary 
conditions for the parameters in Eq. (^) are given by (rn\^)j = (Tn\) 1 - = rn^Sj, ml = 
rrij- = itlq, A 1 ^ = J^Aq, and A % 1 = f^Ao at the Planck scale. Two of the five parameters 
(usually, B and are fixed to obtain the desired values of the Z-boson mass and tan/?. 

In this analysis, we fix the top Yukawa coupling at the Planck scale to be ft(M) = 2.4 
which corresponds to, at the GUT scale, ft(M G ) = 1.4. We choose this value partly for 
comparison with Ref. 0]. It gives the top mass m t ~ 170 (180, 190) GeV for tan/? = 2 
(3, 10). Also we choose A = 0. Dependence of the size of the generation mixing in the 
slepton masses squared on these parameters was discussed earlier. 

Branching ratio for /1 — > e^ is shown in Fig. 3. The horizontal line in the figure is taken 
to be the physical right-handed selectron mass mg R , and M 5 is fixed such that the Bino 
mass Mi is 50 GeV. In our numerical analysis, we imposed theoretical and experimental 
constraints, i.e., negative searches for the SUSY particles at LEP1 and Tevatron |TTJ, a 



- 8 ), 



muon anomalous mag netic dipole moment @ (-1.34 x 10" 8 < \(g - 2) < 2.34 x 10 

and tree-level stability of the scalar potential along h\ = h,2 direction. We present results 

in the case of tan/? = 3, 10, 30 for each sign of \i. 

As can be seen from this figure, there exists a region where the cancellation between the 

diagrams significantly reduces the event rate of fi — > ej. With the universal scalar mass 

hypothesis, an absolute value of Higgsino mass /i determined by the SU{2) L x U(1)y 

^The Yukawa coupling constant we adopt may result in a top quark mass which is fairly large. However, 
we shall note here that the event rate becomes smaller for a smaller value of ft- In that case, experimental 
verification of LFV becomes severer than the case discussed in this paper. 
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symmetry breaking condition increases with in a range 180GeV £ |/x| £ 200GeV 
for tan 0=10 or 30, Mi = 50 GeV and m iR < 300 GeV (200GeV < |/x| < 240GeV for 
tan 0=3). Since the mass range for /x is comparable with other superparticle masses, we 
cannot ignore the diagrams which have the Yukawa type vertices. In a small region 
A\ |b tends to dominate over other contributions, and in a large m gi j region A\ | c does, 
especially if tan/3 is large. Therefore, in a middle, there is steep cancellation. If the 
Yukawa coupling of Higgsino were ignored, the branching ratio would exceed 10~ n in the 
case, for example, when tan (3 = 30 and mg R < 200GeV. Under the universal scalar mass 
hypothesis, however, the cancellation between diagrams reduces the branching ratio to 
~ 1O~ 12 .0 Also, in the region of valleys A\ | a and A\\& contribute to the amplitude with 
the same order as the sum of A\ |b and A\\ C1 and this leads to the difference of positions 
of valleys of each line. Moreover, if the top Yukawa coupling is smaller, the radiative 
symmetry breaking condition requires smaller values of |/x|, and the point of cancellation 
tends to shift to a region with smaller right-handed selectron mass. Thus, the situation 
is more complicated than that previously considered in Ref. |4[|. 

Next, we show branching ratios for /x — > eee, /x-e conversion in ||Ti, and r — > /X7 in 
Figs. 4-6. Configurations of these figures are the same as Fig. 3. These events are reduced 
by about one order of magnitude compared to the case where the Yukawa vertices are not 
taken into account, and these processes also suffer from the partial cancellation. 

The behavior of the /x — > eee is similar to that of /x — > e^ up to normalization. This is 

because the penguin type contribution from Eq. ([0]) is enhanced by phase space integral 

and hence it dominates over other contributions ||10[- Steepness of the valleys is somewhat 

less in the case of /x — > eee, because diagrams other than the penguin diagram contribute. 

On the other hand, the enhancement by the phase space integral does not occur for the 

/x-e conversion and thus behavior of this event rate is not similar to that of /x — > ej. Then, 

§ The event rate of /x — ► cy does not drop off as quickly as the fourth power of the right-handed 
selectron mass above the cancellation point, but has rather a flat behavior. Since the renormalization 
effects by the top and bottom Yukawa couplings reduce the right-handed stau mass significantly, the loop 
diagram like Fig. 2 (c), which is proportional to (m~ )^(m~ )£, dominates the other diagrams completely 
when the right-handed selectron mass is larger than the bino and Higgsino mass. While this fact does 
not change the qualitative behavior that there is a cancellation between diagrams as we mentioned in 
the previous footnote, the quantitative behavior in that region changes even with right-handed selectron 
mass larger jl3|]. The other LFV events also have such smooth behavior as that of fj, — > erf. 
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the cancellation occurs at right-handed selectron mass different from the case of fi — > 
and this suggests search of fi-e conversion would have a complementary role with that of 
jj, — > ej. Also t — ► fi'j behaves like [i — > ej does and they are related with each other 
through the KM matrix. 

So far, we have assumed the universal scalar mass at the Planck scale. It is well-known, 
however, the degeneracy of all soft scalar masses is not needed to avoid disastrously large 
flavor changing neutral currents. In the context of SU(5), the mass degeneracy of the 
sfermions with the same quantum numbers (5* and 10) is sufficient. Thus we shall discuss 
the fi — > e7 rate as an example in a non- universal case. 

To study the case with non-universality, we show how the results depend on the mass 
spectrum of the superparticles with a fixed value of the right-handed mass matrix 



(These values are realized for the universal mass case with m e - K = 300 GeV, Mi = 50 
GeV, A = 0, tan/? = 3 and ft(M) = 2.4.) In Figs. 7 and 8, branching ratio for the 
process [i — > ej is shown as a function of the physical left-handed selectron mass mg L and 
/x-parameter. Here we take Mi = 50 GeV, mg R = 300 GeV, and tan (3 = 3 (Fig. 7) and 30 
(Fig. 8). We assume that all left-handed charged slepton have a common mass mg L and 
A(rriz) = for simplicity. The shaded regions are excluded by the present experiments. 

In general, the amplitude for the fi — > ey is a complicated function of parameters 
at low energy scale. In addition to the right-handed slepton masses squared, it depends 
on the following parameters at the electroweak scale: wij L , A, fi, tan/3, Mi (and M2). 
However, the qualitative behavior can be easily understood by using the mass insertion 
formula discussed above. For example, the left-handed slepton mass becomes smaller 
or Higgsino mass fi larger, A\ |b in Eq. (|T4]) becomes larger than the others. When the 
left-handed slepton mass becomes larger, A\ | a and A\ |d ( A\ | c ) becomes dominant for 
tan j3 = 3 (tan /3 = 30). Thus, the branching ratio becomes larger because the cancellation 
does not occur significantly. 

We should emphasize that our low-energy approach given here should apply to a 
certain class of the unified models where the right-handed slepton masses are the only 



2 2 
m r = m- 



f 1.0 -9.8 x 10- 5 - 

-9.8 x 10~ 5 1.0 
\ -3.2 x 10~ 3 -2.3 x 10~ 2 



3.2 x 10~ 3 

2.3 x 10~ 2 
0.25 




(19) 
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source of flavor violation. The effect of this interference can be important when one 
computes the event rate in a given model and we expect that the rates suffer from this 
complicated cancellation in a large class of models. On the other hand, there is not 
such significant cancellation in models which left-handed sleptons have LFV masses, e.g., 
SO (10) SUSY GUT 1, |, | and the minimal SUSY standard model with right-handed 
neutrinos [0, 0, ITOl since the LFV process is usually dominated by only one diagram. 
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Figure 1: Top quark mass dependence of flavor- violating right handed slepton mass (to- 
at the electroweak scale. Real line, dashed line, and dash-dotted line correspond to the 
case for tan/5 = 30, 10,3. Here we take M 1 = 50 GeV and TOg fl = 300 GeV, A = 0. 
We represent the points corresponding to ft{M) = 0.5,1.0,2.0,2.4,5.0, which is the top 
Yukawa coupling constant at the reduced Planck scale. 

Figure 2: Feynman diagrams contributing to the process /x — * &y in the minimal SUSY 
£77(5) GUT. The symbols JIr, Br, B°, and H° represent left-handed smuon, right- 
handed smuon, right-handed selectron, Bino, and neutral Higgsino respectively. The 
blobs in the slepton line indicate the insertion of the flavor-violating right-handed slepton 
mass (to- ){f and left-right mixing mass which is proportional to the vacuum expectation 
value (hi) or (h 2 ). The blobs in the Bino-Higgsino line represent the insertion of gaugino- 
Higgsino mass mixing, that is, \i denotes Higgsino (Hi — H 2 ) mass mixing, (h^ 2 ) Higgsino- 
Bino (H° — B°) mass mixing, Mi Bino mass. And blob in the external muon line represent 
the chirality flip of the external muon /i. 

Figure 3: Branching ratios for the process fi — > cy as a function of the physical right- 
handed selectron mass, mi R . Real line, dotted line, dash-dotted line correspond to the 
case for the tan j3 = 30, 10, 3 respectively. Dashed line represents the present experimental 
upper bound for this process. Here we take Mi = 50 GeV, ft(M) = 2.4, and (a) /j, > 
(b) n < 0. 

Figure 4: Branching ratios for the process /i — > eee. Here we take the same parameters 
as Fig.3. 

Figure 5: The /i-e conversion rates in nuclei f^Ti- Here we take the same parameters as 
Fig.3. 

Figure 6: Branching ratios for the process r — > /ry. Here we take the same parameters as 
Fig.3. 

Figure 7: Branching ratios for the process \x — > as a function of the physical left-handed 
selectron mass mg L and /i-parameter for the case of non-universal scalar mass. Here we 
take Mi = 50 GeV, m~ eR = 300 GeV, f t (M) = 2.4, tan (5 = 3, A(m z ) = and the typical 
value of right-handed slepton mass matrix as mentioned in text. For simplify we assume 
that all left-handed charged slepton have a common mass rrie L . The shaded regions are 
excluded by the present experiments. 

Figure 8: Same as Fig. 7 except for tan/5 = 30. 
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